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Abstract 


In this note we show that the 
Fourier approximation to the operator 


eigenvalues 

sin(2x) 


of the pseudospectral 
satisfy 


R e z i = ± 1 or R e Z i = °- 


Whereas this does not prove stability for the Fourier method, applied to the 
hyperbolic equation 

U t = sin(2x)U x - ir < x < u; 

it indicates that the growth in time of the numerical solution is essentially 
the same as that of the solution to the differential equation. 
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1. Introduction 


Let us consider the problem 

U - GU = 0 0 £ x £ 2tt 

U(x,0) = U°(x) (1.1) 

where 

G = a(x)|^ . (1.2) 

In the Fourier pseudospectral (collocation) method, we seek a trigonometric 
polynomial of degree N , U^that satisfies 

- g n u n - 0 

U N (X,0) = U°M (1.3) 

where 

g n = p n g ; 

is the pseudospectral projection operator [5]. It is known [2] that 
when a(x) does not change sign in the interval, the semidiscrete solution 
of (1.3) is stable. When a(x) changes sign in the interval, the 
situation is much more complicated. Gottlieb, Orszag and Turkel [1] have 
proved stability for the case where a(x) is of the form 

a(x) = a sin(x) + 6 cos(x) + y . (1.4) 

In [4] , Tadmor argues that this stability proof results from the special 
form of a(x) in (1.4) and cannot be extended. In the next section we 
prove a theorem related to the problem of stability of (1.1) where a(x) 
is a second degree trigonometric polynomial. 
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2. The Theorem and Its Proof 

Theorem : Considering (1.1), (1.2), where a(x) = sin(2x), then the eigen- 

N 

values X. of G x , satisfy 
1 N 

Vi - - 1 or R e x i = 0 or Vi = 1 - (2.1) 

Proof : 

The projected subspace that results from using the operator 

is spanned by the following 2N basis functions 

V N = ^{1 ,cos (x) , . . . , cos (Nx) , sin(x) , . . . ,sin(N-l)x) } t (N even) (2.2) 

Define the following four subspaces of V 

N 

Wj = S p {cos(x), cos (3x) , . . . ,cos((N-l)x) } 

W 2 = Sp{sin(x) , sin(3x) , . . . ,sin((N-l)x) } (2.3) 

W_ = S (sin(2x, sin(4x) , . . . ,sin((N-2)x) } 

3 P 

W 4 = S p^ 1 > cos ( 2x )>-*-» cos ( Nx )). 

It is easily verified that 

V N = W 1 ® W 2 ® W 3 ® W 4 (2 .4) 

and each is invariant of G^; therefore we can discuss separately 

the four matrices which represent G N in each one of the subspaces HL. 
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Define now 

B M = [G„l 1 $ i $ 4 (M=|); (2.5) 

l L N £ 

M 

then by using elementary trigonometric relations we get that are 

tridiagonal matrices whose elements are: 





-5- 


let A by any tridiagonal matrix: 


b 2 a 2 


b . a c 

n-1 n-1 n-1 

b a 

n n 


(2-6) 


and let by the submatrix 


\ " 


b k-l a k-l C k-1 
b k a k 


(2.7) 


Upon defining 


q k (A) - det A r 


( 2 . 8 ) 


it is easily verified that 


W A) ■ Vi <1 k (A) - b k+i c k ’k-i (A) 


(2.9) 
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and 

q R (A) = det A. 

In the following we treat each one of the matrices B^, i = 1,2, 3,4 
separately. 

Lemma 1 : The matrix B^ has one zero eigenvalue, and all its other 

eigenvalues X i satisfy R g X^ = 1. 

Proof : For any M define 

Si * 2B i ' XI • 

The characteristic polynomial of 2B^ is given by 

^(X) = det Si (2.10) 

and using (2.8) 

V X > - W- 

We define now the following family of polynomials (in the variable X) 


P 0 = 1 Pj = - (X + 1) 

p k.i ■ - xp k * < 4k2 - U P k-1 1 5 * * * • 

Note that from (2.9) and the structure of C. 


( 2 . 11 ) 


P k = W 2 s k < M ; 


( 2 . 12 ) 
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* 

however (2.12) is not true for k = M; rather we have 

Q m 00 = (2M - 1 - X) P M _ 2 + (4 (M-l) 2 - 1) P M _ 2 2 < M . 

From (2.11) we get 

' (2M - 1) P M1 + P M 2 < H. 

Using (2.14) and (2.13) results in 

Vl M J P M * 2 < M - 

Finally we solve (2.15) for P M in terms of Q^fA), Q^ +1 (A) and 
substitute the result in (2.14). We thus get the polynomials Q^(A), M 
that satisfy the following recursion formula 

Q 2 0) = X(X-2) ; Q 3 (X) = - X(X 2 -4X + 13) 

Qm + i(U = (2-X) <^(X) + (2M-1) 2 Q m _ 1 (X) 3 < M . 

It is easy to verify now that X = 0 is an eigenvalue of 2B' 1 . In fact 
X = 0 is a root of Q 2 (X) and Q 3 (X) and therefore of any Q M (X). We 
define now 

x = i(2 - X) (a) 

and 

■VjW - r QmW • w 

to get 

R 2 = - x ; R 3 = x 2 - 9 

and 

Vl ■ * - C ^- 1 ) 2 “M-l 


(2.13) 


(2.14) 


(2.15) 


2 


(2.16) 


(2.17) 


M :> 3 . 


(2.18) 
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The relation (2.18) defines R^fx) as a family of orthogonal polynomials 

on the real axis. Therefore, for every M the roots of R^x) are real, 

which implies by (2.17) (a) that 2 - X are imaginary. Therefore, the 

M 

eigenvalues of the matrices 2B^ for any M have real part equal to 2. 
This completes the proof of Lemma 1. 

M 

Lemma 2 : For any M the matrix has one zero eigenvalue and the real 

part of the others is -1. 

Proof : The proof is an immediate result of the fact that in view of (2.9) 

q k C-B^ - xi) 

M 

satisfy the same recurrence formula as 

M 

Lemma 3 : The eigenvalues of are purely imaginary. 

Proof: Define the matrix 

1//2 

1//4 

lM^2 

Then it is clear that 

D ' 1 Bj D 

is a skew symmetric matrix, and therefore its eigenvalues are purely 

- M 

imaginary. The same is of course true for B^. 

M 

Lemma 4: The eigenvalues of B 4 are purely imaginary. 
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M M 

Proof: From the definition of B„ and B. it follows that if P. is 

3 4 k 

M 2 

characteristic polynomial of then X is the characteristic 

polynomial of (B 1 ^) ^ + 2) x (k+2) * T ^ us t * ie ei £ envalues °f B 4 are purely 
imaginary. 


The proof of Lemma 4 concludes the proof of the theorem. 
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